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We find a family of exact solutions of Einstein equations describing the field of a static axisym- 
metric thin disk living in six-dimensional space-time. In particular, we study the disks constructed 
(by cutting out the central part of the space-time) from the conventional Schwarzschild and Chazy- 
Curzon solutions with simple extensions in the extra 2 dimensions. The disks are interpreted in 
terms of two counter-rotating streams of particles on free circular orbits. Two extra parameters - 
the constants of motion resulting from projections of the particle-velocities' extra components - are 
constrained by the requirement that the orbits within the disc be stable. The requirement is met 
just in cases when the radial profile of the disc orbital speed fits the rotation curves of many spiral 
galaxies. It thus turns out that the effective modification of the gravity law by extra dimensions 
could explain the observed flatness of these curves equally well as the usually given dark matter 
interpretation. In the second part of the paper we show that the inclusion of extra dimensions also 
leads to better fits of the gravitational lensing data for galaxy clusters, without changing results 
obtained in solar system scales. Finally, we discuss whether the effect of extra dimensions could 
also be translated as following from the occurrence of extra matter ("Kaluza-Klein particles"). A 
comment on possible properties of such particles and a chance to detect them at LHC is given in 
the appendix. 
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INTRODUCTION 



In the present work we are interested to model a grav- 
itational 4D disk living in a universe endowed with ex- 
tra dimensions. The motivation comes from the various 
properties one can derive as rotation curves and density 
and pressure profiles now carried on with a more general 
geometrical approach for space-time. A first serious ac- 
complishment at this direction is only achieved by investi- 
gating solutions of Einstein field equations in axially sym- 
metric configurations in D dimensions. Exact solutions 
of Einstein field equations in axially symmetric configura- 
tions are an important tool to understand the dynamical 
properties of real systems which can be described approx- 
imately by a thin disk. In this sense, axially symmetric 
solutions have great astrophysical interest, because they 
can be used to model galaxies and accretion disks. A 
long range of disk solutions was derived with or without 
radial pressure. Solutions for static disks without radial 
pressure were first studied by Bonnor and Sackfield 
and Morgan and Morgan [2|, and with radial pressure 
by Morgan and Morgan A large class of static thin 



disks solutions were obtained by Letelier and Oliveira 
Q using the inverse scattering method. Disks with ra- 
dial tension have been considered in [H], and disk models 
with electric fields 0| and magnetic fields and both 
magnetic and electric fields @]. Solutions for self-similar 
static disks were analyzed by Lynden-Bell and Pineault 
01, and Lemos [13] • Another relevant approaches are 
the superposition of static disks with black holes or disk 
solutions with halos [ll|, and see other important astro- 
physical solutions to mimic AGNs in Recently Vogt 
and Letelier refined general relativist models of galaxies 
in a considerable way [l3j]. Important discussions about 
the role of general relativistic disks to explain rotation 
curves of galaxies can be found, e.g., in [14J and essential 
counter-arguments in [ljj]. On axisymmetric geometries 
in D dimensions, generalized Weyl solutions in D dimen- 
sions were presented by Emparan and Reall [16| and im- 
portant objects were classified by them, as black rings 
and Kaluza-Klein bubbles. About disk rotation curves 
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in higher dimensions see for instance [17j. For other im- 
portant works about solutions of Einstein equations in D 
dimensions see (l8| . 

The treatment proposed at the present work permits 
to find important imprints of extra dimensions in the 
density profile of the disk. A first result obtained by 
introducing extra dimensions shows important modifica- 
tions in rotation curves profiles. As rotation curves are 
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the major tool for determining the distribution of mass 
in disk galaxies, and are also important to study kine- 
matics and to infer the evolutionary histories in galactic 
systems, as also the most basic and classic manner to 
bring on the presence of dark matter in galaxies, maybe 
it is possible that a relativistic disk model embedded in 
a multidimensional universe shed some light on the dark 
matter problem. Actually, dark matter (DM) is one of 
the most intringuing problems in astronomy, astrophysics 
and cosmology. This puzzle has paramount consequences 
on the deep comprehension of the universe, matter in gen- 
eral and particle physics in particular. Rotation curves 
of galaxies and gravitational lensing of galactic clusters 
are the main basic observational tools to understand the 
presence of DM in universe On this behalf, the na- 
ture of DM can be investigated by applying a range of 
gravitational models: by postulating a halo of exotic mat- 
ter around galaxies (e.g. CDM models [2(3]) or by mod- 
ifying grav ity at large scales (e.g. MOND, TeVeS and 
STVG [21 1; a s an important review about the subject 
see, e.g., 22J). Another relevant alternative approach 
is to enhance astrophysical and cosmological models by 
using extra dimensions 23, 24, 2^, If]]. Here, to under- 
stand the role of dynamics in producing a DM behavior in 
galactic disks, as a first step, we model and investigate 
a general relativistic disk and also calculate the gravi- 
tational lensing regime for a galactic cluster living in a 
universe endowed with extra dimensions. We adhere to 
the standard general relativity, just in space-time with 5 
spatial dimensions. 

An important question, unanswered until now, is why 
a galactic disk does not follow a gravitational Newto- 
nian or general relativistic 4D profile for rotation curves. 
The extradimensional argument is a proposed solution in 
particle physics field: there the cosmological amount of 
dark matter is used to calculate cross sections of postu- 
lated Kaluza-Klein (KK) particles to be tested in future 
experiments at colliders [23]. Here we follow another di- 
rection to show quantitatively how one can calculate ro- 
tation curves from GR in a higher dimensional pattern. 
Therefore, we construct a disk endowed with extra di- 
mensions by solving the vacuum Einstein equations for 
an extension of the Weyl's metric. The main purpose of 
the work is to construct a top-down model, where we are 
interested in what is coming to pass at large scales as- 
suming a multidimensional scenario. Usually thin disks 
are solved from Chazy-Curzon or Schwarzschild solutions 
and as the simplest example, we solve the case for a 6D 
space-time. The system is constructed from such com- 
monly used solutions: they are employed to generate the 
disk by introducing discontinuities in the first derivatives 
of the metric functions (image method) . Two constants of 
motion from projection of extradimensional particle ve- 



locities are the free parameters of the model. We prevent 
the ad hoc adjustment of such parameters with observed 
rotation curves, preferring to investigate values where the 
disk becomes stable. The stability is achieved when the 
disk is Newtonian-like (where such parameters are null) 
or for a tiny range of values that astonishingly makes 
the circular geodesies fit with great precision the rota- 
tion curves of many spiral galaxies. We compare such 
results to well succeeded astrophysical dark matter pro- 
files and demonstrate that our predictions recover natu- 
rally gravity at solar system and give the same gravita- 
tional lensing as does a dynamically successful dark halo 
model. We propose, following such results, that the ge- 
ometry of a space-time endowed with extra dimensions 
modify gravity at certain distance/mass scales and dark 
matter is replaced by the picture of a multidimensional 
universe. We also consider, for completeness, another 
interpretation where the geometry can induce Kaluza- 
Klein fields and particles on the AD space-time, and the 
present model could be a interesting constraint for theo- 
ries on KK dark matter particles. 

The present work is organized as follows: in Section ITT1 
we present the principal aspects of static thin disk solu- 
tions in AD and the expression for the energy-momentum 
tensor of the disk. In Section IIII Al we derive general 
Weyl's solutions. To derive Weyl's solutions to higher 
dimensions is to find all solutions of the vacuum Ein- 
stein equations that admit D — 2 orthogonal commuting 
Killing vector fields. We show, following [if]], that in 
D dimensions, Einstein equations in a Weyl geometry 
is reduced to the Laplace's equations, what permits to 
use the image method to construct a disk living in a D 
dimensional universe. In Section IIII B I we exemplify by 
constructing a simple multidimensional disk model using 
six dimensional axially symmetric Einstein equations in 
vacuum, with a "Schwarzschild solution" for the four di- 
mensional part; the extra dimensional part is solved with 
a Chazy-Curzon solution. These solutions are employed 
to generate a disk by introducing discontinuities in the 
first derivatives of the metric functions (image method) . 
In Section llVl we find a equation for the circular geodesic 
orbits (rotation curves) and study the stability of such or- 
bits in Section [V] both by the Rayleigh criterion and by a 
perturbative method. Results and graphics for stable ro- 
tation curves and comparison to observed rotation curves 
of spiral galaxies are showed in Section IV CI In Section 
IVII we introduce quantitative arguments to show that at 
solar system scales, the model presents the correct be- 
havior for rotation curves (a Newtonian-like profile). In 
Section IVTT1 we derive the gravitational lensing regime of 
the model to a spherically symmetric cluster of galax- 
ies. Discussions about such results is done in the Con- 
cluding Remarks (Section I VIIip . And although the main 



purpose of the present work is to construct a top-down 
model and not to shed light on a fundamental microscopic 
theory we also illustrate another possibility, considering 
that our model could constrain a Kaluza-Klein dark mat- 
ter particle to be tested at Large Hadron Collider (LHC) 
in next years, where we write considerations about the 
compactification of a 6D space-time following the exam- 
ple presented. In what follows, we do c = G = 1. 



II. THIN DISKS IN 4D 

Much endeavor has been strongly attached to advance 
in techniques for finding gen eral relativistic exact solu- 





FIG. 1: Illustration of the "displace, cut and reflect" method 
for the generation of disks. In A the spacetime with a singu- 
larity is displaced and cut by a plane (dotted line), in B the 
part with singularities is disregarded and the upper part is 
reflected on the plane. 



tions in four dimensions [2jj, \M- One of the earliest Noting that W(Q = W((), one can write dW^dW^ 



accomplishments in this direction was attained by Weyl 
[29! ]. who found the general static axisymmetric solution 
of the vacuum Einstein equations: 



ds 2 



e 2 *dt 2 + e- 2 * (e 2 ^(dr 2 + dz 2 ) + r 2 d/) , (2.1) 



where <£>(r, z) is an arbitrary axisymmetric solution of 
Laplace's equation in a three-dimensional flat space with 
line element 



da 2 = dr 2 



2 d V 2 



dz 2 , 



and 7 satisfies 



7,r 



7,; 



2r$ r $ 



(2.2) 



(2.3) 



(2.4) 



where ( ) jQ = d/dx a . The solution of these equations 
is given by a line integral. Since $ is harmonic, it can 
be regarded as a Newtonian potential produced by cer- 
tain (axisymmetric) sources. Since in this coordinate the 
spherically symmetric solutions of the Einstein equations 
correspond to a bar of density 1 /2. one needs to be careful 
in the use of Newtonian images [4|. 

As an illustration, the metric l|2.ip for a static axially 
symmetric AD space-time can be written in quasicylin- 
drical coordinates (r, <p, z) in the form 



ds 2 



-'''dt 2 + X 2 e 4 'd<p 2 + f(dr 2 + dz 2 ), (2.5) 



where x, 4>, and / are functions of r and z only. In the 
vacuum, the Einstein equations are equivalent to 



X,rr + X,zz = 0, 



(X0,r),r + (x<f>,z),z = °- 



(2.6) 



(2.7) 



^-^dCdC = |VF(C)| 2 dCdC. Or even, dWdW = d X 2 + 
dZ 2 = |W(C)| 2 (dr 2 + d;z 2 ). Thus without lossing gener- 
ality one can choose x = r - In such a way we can write 

f(r, z) as 

2 . 

(2i 



ln/[0] = 



m- 

[24>^. z ]dz}. 



h 2 z ]dr 



In order to obtain a solution of i|2.6p - i|2.8p representing 
a thin disk located at z = 0, we assume that the met- 
ric functions x, 4>, an d / are continuous across the disk, 
but have discontinuous first derivatives in the direction 
normal to the disk. Thin disk solutions in Weyl coor- 
dinates are functions of the class C° . The reflectional 
symmetry of l|2.6p - (|2,8|) with respect to the plane z = 
allows us to assume that x, tfi, and / are even functions 
of z. Hence, x,z, 4>,z and f ;Z are odd functions of z. We 
shall require that they not vanish on the surfaces z = . 
Such impositions can induce e.g. a discontinuity in the 
space-time by reflecting the solution through the plane. 
This represents the construction of the disk using the well 
known "displace, cut and reflect" method that was used 
by Kuzmin [3jJ|in Newtonian gravity and later in GR by 
many authors ]H|-[i3|- The material content of the disk 
will be described by functions that are distributions with 
support on the disk. The method can be divided in the 
following steps that are illustrated in Fig. [TJ First, in a 
space wherein we have a compact source of gravitational 
field, we choose a surface (in our case, the plane z = 0) 
that divides the space in two pieces: one with no singu- 
larities or sources and the other with the sources. Then 
we disregard the part of the space with singularities and 
use the surface to make an inversion of the non-singular 
part of the space. This results in a space with a singu- 
larity that is a delta function with support on z = 0. 



Let ( = r + iz. It is possible to consider \ as the real 
part of an analytical function W(() = x(r, z) + iZ(r, z). 



This procedure is mathematically equivalent to make 
the transformation z — > \z\ + a, with a constant. In the 
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Einstein tensor we have first and second derivatives of z. 
Since d z \z\ = 2d{z) - 1 and d zz \z\ = 2S(z), where <&(z) 
and 8{z) are, respectively, the Heaviside and Dirac distri- 
butions, the Einstein field equations will separate in two 
different pieces [3J|: one valid for z ^ (the usual Ein- 
stein equations), and other involving distributions with 
an associated energy-momentum tensor. Due to the dis- 
continuous behavior of the derivatives of the metric ten- 
sor across the disk, the Riemann curvature tensor con- 
tains Dirac delta functions. The energy-momentum ten- 
sor can be obtained by the distributional approach due 
to Papapetrou and Hamouni [j^ ]. Lichnerowicz jH3|, and 
Taub [3. It can be written as T a b = [T a b ] S(z), where S 
is the Dirac function with support on the disk and [T a b ] is 
the distributional energy-momentum tensor, which yield 
the volume energy density and the principal stresses. A 
different approach is given by Israel [35J where one makes 
use of the extrinsic curvature of the surface to represent 
the matter. Here, a, b = 0, 3. 



A. Density and pressures on the disk 

The disk at z = divides the space-time into two 
halves. The normal to the disk can be described by the 
co- vector n a = dz/dx a = (0, 0, 0, 1). Above the disk near 
z = 0, we can expand the metric as 



n t — n° 4- ~ I _i_ „2 ^jab I 

9ab - 9ab + Z-7^-U=0 + z -q^t\ 



and below z = 0, 

9ab = glb + z 



dz 



dz 2 



(2.9) 



(2.10) 



The quantity g® b means the value of g a b at z = 0. The 
discontinuities in the first derivatives of the metric tensor 
can be cast as 



9ab,; 



9ab,. 



(2.11) 



The symmetry of the problem gives 4> + z \ z=0 = -^|, =0 , 
and Z+Uo = -f~\ z=0 . Denoting 6. z \ z=Q = (t> + z \ z=0 , and 
the same for the other functions, we calculate the follow- 
ing discontinuities: 



2e~ 



brr = 2f. 



(2.12) 



(2.13) 



It is possible to find the other discontinuities terms by 
contracting indices in b ab : 



Z|z=0 ' 



b% = 26 z , 



h rr = — — f I 



p2 $ > z I s=o 



b VV 



2 f f,z, 



~f,Z ) 



-20,; 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



From (|2.1ip . one can work out (the discontinuities of) 
the Christoffel symbols through the disk given by 



[r a 6 



-(b a c 5 z b + b a b 6 z c -g az b bc ) 



(2.20) 



where [r a 6c ] = r+ a fcc - T- a bc at z = 0. From the Rie- 
mann tensor defined by 

Rabcd — ~^{.9ad,bc 9bd.ac 9bc,ad 9ac,bd) 

+9rsF r adF S bc ~ 9rs^ r ac^" bd, (2-21) 

we can compute the Riemann distributional tensor, 

[R a bcd] = \{b a d 8\5 z c - b\5 z b 5 z d + g az b bd S z c ). (2.22) 

Defining the Ricci distributional tensor as [R a b] = [R c acb] 
and the Ricci distributional scalar [R] = [R a a ], we can 
identify the distributional energy-momentum tensor on 
the disk through Einstein equations as 

[R a b ] - ^6 a b [R] = 8tt[T\}. (2.23) 

Then the distributional energy-momentum tensor is 
given by 

\rjia i t uaz cz is2 m , az-iz zzia 

[T b\ = J£-{b o b-b d b +g b b -g b b 



+b c c (g zz 6 c b -g az 6 z b )}. 



(2.24) 



This energy-momentum tensor describes the matter con- 
tent (fluid) of a thin disk located on z = 0. The compo- 
nents of a such tensor are calculated as 



1 



[T t t} = —{-b* z +g**(b r r + b* z )}, 



(2.25) 



&zz — hrr-. 



(2.14) 



lT r r] = T ^{-b zz +g zz (b t t + b%)}, 



(2.26) 



2rV0 ; 



(2.15) 



[T z 



(2.27) 
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16tt 

Defining the vierbein 



(2.28) 



e [t) a = (-^=,0, 0, o) , e {r) a = ( 0, -L, 0, o) , 
e (v) a = (0,0,— ^,oV e (2) a - [0, 0,0,- 



one can write down the energy-momentum tensor <|2.24|) 



a b , a b 

ee(t) e (t) + p r e( r ) e (r) 



+ Pv e (v) ae (v) & +Pze( z ) a e (z) , (2.29) 



a „ 6 



yielding the volume densities, i.e., the energy density and 
pressures as 



= -[T\] = 



2 l»=o 



P, = [T%]=-^\ 



Pr = [T r r ] = 0, 



Pz = [T z z ] = 0. 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



As an example, in Fig. [2] we show the density e and 
pressure p v profiles of the AD constructed thin disk for 
a Schwarzschild solution (as explained bellow), see also 



B. Solutions 



C. Rotation curves 

In first approximation one can consider that the par- 
ticles of such above fluid move along geodesies. In par- 
ticular, we can consider particles moving along circular 
geodesies whose tangential velocities give us the rotation 
curves. 

From Eq. I|2.5p we have the first integral of motion, 

-e-^ 2 + /(r 2 + i 2 )+r 2 cV = 1, (2.35) 

where x a = dx a /ds. Assuming f = and z = (particles 
with no radial motion and confined on z = 0), Eq. (|2.35p 
reads 



-All i 2 <b -2 -i 

e v t + r e v Lp = 1. 



(2.36) 



The geodesic equations on the disk reduce to 

(e-^i 2 - (rV). r (^ 2 = 0. (2.37) 

Eqs. (|2.36j) and (|2.37p form a system of equations for Cp 1 
and t 2 . From these equations we find the rotation curves 
Vc, 



_ 9 ^ = _ 9 ^ 
9tt at V 9tt t 2 



reduce to 



V c = v /-0, J ./(2/r + 0, r ). 



(2.39) 



In Fig. [3] we show some rotation curves for the con- 
structed AD disk, where we are varying the cut disk pa- 
rameter a. Here we are applying the thin rod solution 
with linear mass density of 1/2 (Schwarzschild solution). 



In Eq. (|2.5p . the function 4> ls intrinsically related to 
the Newtonian potential U by 4> — 2Z7. A important 
property of Weyl's metric is the fact that Eq. I|2.7p to 
be a Laplace's equation in cylindrical coordinates, and 
exploiting the characteristics of its linearity it is possible 
to employ the superposition of solutions. Some of the 
most commonly used asymptotically plane solutions are, 
e.g., Chazy-Curzon, the finite rod, and the Schwarzschild 
ones. 

In the Chazy-Curzon case, the solution for_a particle 
of mass m in the position z = z is given by 



36 



2 9 

m r 



2m 

~R' ln/= B* 



(2.34) 



where R = ^Jr 2 + (z — z ) 2 . 

The Schwarzschild solution corresponds to taking the 
source for to be a thin rod on the z-axis with 1/2 linear 
mass density. 



III. AXIALLY SYMMETRIC SOLUTIONS FOR 
ANY D AND THE DISK CONSTRUCTION 

A. Weyl's geometry for D dimensions 

The first action to generalize Weyl's geometry to space- 
times endowed with higher dimensions is to consider a 
suitable coordinate chart for the D-dimensional line el- 
ement. Here we follow the method presented by [lfjj], a 
simple generalization of what is done in four dimensions 
[2^ |. It will be alleged that the metric is Riemannian or 
Lorentzian. Let £u) represent the Killing vector fields, 
I < i < D — 2. It is admissible to elect coordinates 
{x' l 1 y 1 1 y 2 ) such that = d/dx l , since these commute, 
with the metric coefficients depending only on y 1 and y 2 . 

Now we must show that one can select coordinates y 1 
and y 2 to span two-dimensional surfaces orthogonal to all 
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FIG. 2: (a) AD disk surface density per unit mass profiles for disk cut parameters a — 1, 1.5, 2, 2.5, 3 (from top to bottom). We 
take r' = r/m. (b) The same as (a) but for disk pressures. The solution used here is the Schwarzschild one (Chazy-Curzon 
also offers similar profiles). 
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where Z and Z are complex conjugate coordinates if the 
transverse space is spacelike. The function C is inde- 
pendent of x % . Now, assuming that \i,v = i,Z,Z, then 
the vacuum Einstein equations read = 0. The ij 
component gives 



exp [ d z^ 



exp dz$z 



0. 



(3.42) 



Summing this equation over i yields 



FIG. 3: AD disk rotation curves with cut parameters from 
a — 0.3 to a = 2 (from top to bottom). These curves have a 
Newtonian-like profile. We take r' = r/m. The solution used 
is the Schwarzschild one (Chazy-Curson also offers similar 
profiles). 



of the . In order to perform this, one has to make ev- 
ident that the two-dimensional subspaces of the tangent 
space orthogonal to all of the vectors are integrable. 
Sufficient conditions for integrability are afforded by the 
Emparan-Reall theorem [if]]. If it is further assumed that 
the Killing vector fields are orthogonal to each other then 
the metric must take the form 



D-2 



(3.40) 



where a and b take the values 1,2, the metric coefficients 
are independent of x z , and = ±1 according to whether 
is spacelike or timelike. Locally it is always possible 
to find coordinates such that 



gabdy a dy b 



,2C 



dZdZ, 



(3.41) 



(3.43) 



This last has as general solution 

■ = log (w(Z)+w(Z)), (3.44) 

3 

where w = w if Z and Z are complex conjugate. Substi- 
tuting equation (|3,44|) into equation (|3.42|) yields 

2 (to + w)d z d z <$>i + dzwdz^t + d z wd z ^ l = 0. (3.45) 

If to is non-constant then Rzz = can be rearranged to 
give 



Ei<j dz^dz^ 



An akin equation is obtained from R zz = 0: 



(3.46) 



(3.47) 
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Integrating the first two terms of these equations and 
applying equation (|3.44p we get 



c = 2 lo S (dzwdzw) + 3, 



(3.48) 



where 



dz5 = - 



w + k; 



^Sz^iSz^, (3.49) 



%3 



W + W sr-^. „ , _ , 



Kj 



The integrability condition for 3 is 
d z d z E = d z d z E. 



(3.50) 



(3.51) 



It is straightforward to check that this equation is indeed 
satisfied if equations (|3.44|) and l|3.45p hold. These equa- 
tions also confirm that the remaining Einstein equation 
R-zz — i s satisfied. 

Since w and w have been assumed non-constant, one 
can perform a coordinate transformation from Z and Z 
to w(Z) and w(Z) in the similar way as in the 4D case 
mentioned at the beginning of Sec. [TTJ In 4D } they are re- 
ferred as "Weyl's canonical coordinates" This yields 



ds 2 



E 



CjC 



2i>, 



(da 



* 2 + e 2 3 dwdw. 



(3.52) 



This coordinate transformation is conformal. Eqs. 
(|3.45|) , l|3.49p and (|3.50|l are conformally invariant so the 
transformation just replaces d z by d = d w and d z by 
8 = dw- Then the solution is determined by the follow- 
ing equations 



$j = log(w + u5), 



(3.53) 



2(tu + w)<9<9$ 4 + <9$, + 3^ = 0, (3.54) 



Be = -(w + w) d$id<f>j 

i<3 



(3.55) 



(3.56) 



If Z and Z are complex conjugate coordinates then, as 
mentioned above, one must take w — w. Introduce real 
coordinates (r, z) by w = r + iz 7 so the canonical form of 
the metric is 



ds 



e 2H (dr 2 + dz 2 ). (3.57) 



Eq. (|3.54|) then takes the form 
r 9r 



9r 2 



9.: 



0, 



(3.58) 



which is again just Laplace's equation in three- 
dimensional flat space with metric (|2.2|) . The function 
<I>i is independent of the coordinate 9, i.e., it is axisym- 
metric. The solution is therefore specified by D — 3 in- 
dependent axisymmetric solutions of Laplace's equation 
in three-dimensional flat space. 

B. The simplest example: a 6D disk 

The fact that $i, for D dimensions, has as solution ax- 
isymmetric solutions of Laplace's equation leads to the 
possibility to construct a disk with higher dimensional 
components. To obtain a solution for such a metric which 
represents a thin disk located on y a — z — 0, we assume 
the functions of the metric $i and g a b are continuous 
along the disk, in particular on the surface z = 0, but 
with discontinuous first derivatives on that surface. Fol- 
lowing the same approach developed in Sec. HU we can 
introduce these discontinuities by doing the replacement 
z — > \z\ + a, where a is a constant. And as in Sec. [HI 
such discontinuities indicate that the Riemann curvature 
tensor contains Dirac delta functions that makes possi- 
ble to calculate the energy-momentum tensor by a dis- 
tributional approach (see also Sec. [ill for a full expla- 
nation about the method). Calculating the extradimen- 
sional components of the energy-momentum tensor, we 



find that [T£] = -\T$tl] - jk9 zz K>- If D 18 even > 
then the sum of extradimensional pressures cancel and 
the total pressure depends only of AD components. Thus 
we assume a disk living in universe endowed with an even 
number of extra dimensions. Other important arguments 
can be found in literature in favor of an even D (e.g. as 
has been emphasized by different authors, the Huygens 
principle does not hold for odd D [i3|). 

Whatsoever, as the simplest example, we can work out 
in what follows the case for a 6D disk. The metric for 
an axially symmetric 6D space-time can be written in 
quasi-cylindrical coordinates as 



ds 



Vd<r 2 



^+Vc' y dx 2 +c" !/ dy 2 +/(dr 2 +dz 2 ), 

(3.59) 

where (j> = (f>(r,z), f = f(r,z), x = x(r,z), ip = ip(r,z) 
and x and y are the extra dimensional coordinates. The 
vacuum Einstein equations Rab = 0, (A, B = 0, 1, 5) 
reduce to 

(3.60) 



(XV 4'),rr + (xV fp),zz = 0, 



0, (3.61) 
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<l>,r(xVW),r , <t>,z(xVW),z 



x\fi> 



x\fi> 



0. (3.62) 



where ( ), a = d/dx a . About the coordinate free- 
dom, one can consider c = XV^ as the real part of 
an analytical function W(C) = <r(r, z) + iZ(r,z), where 
C = r + iz. Noting that W(Q = W(Q, one can 
write dWdW = ^fpCdC = | V^'(C) | 2 dCdC- Or even, 
dWdW = dc 2 + dZ 2 = \W'(()\ 2 (dr 2 +dz 2 ). Thus with- 
out lossing generality one can choose xVtf = r - We do 
i/j = 1 and x = r i that yields the Einstein equation for 

f(r,z), 

\nf[<b,v] = \Jr{[<j? r -<j> 2 z +24> tr /r + v* r -v%]dr 
+ [2(j> tr( /) tZ +2<f> iZ /r + 2v r v z ]dz}. (3.63) 



Eqs. p.60p - p.63p form the complete set of vacuum Ein- 
stein equations for the metric (|3-59fl . 

To obtain a solution of p.60p - (|3.63p which represents 
a thin disk located on z = 0, we assume the functions of 
the metric / and <f> are continuous along the disk, in par- 
ticular on the surface z = 0, but with discontinuous first 
derivatives on that surface. We introduce these disconti- 
nuities by doing the replacement z — > \z\ + a, where a is 
a constant. The distributional energy- momentum tensor 
reads 



t] = -7-{-b zz +g zz (V r + b z z + b* ip + b x x + ¥ „)}, 

107T 



(3.64) 



1 



r] = T7-{- hZZ + 9 ZZ ( b t + b z z + b* v + b x x + W „)}, 

107T 

(3.65) 



[T z ,] = 0, 



(3.66) 



— {-b zz + g zz (b t t +b r r + b z z + b x x + 6%)}, 

(3.67) 



[t x x ] = _{-&** +<7**(i,* t + b r r + b z z +b* v + v> „)}, 

(3.68) 



t + V r + V , + + x )}. 

167T 

(3.69) 

And the discontinuities of first derivatives of the metric 

b A B = {g\ B .z -9ab.z)\,=o yields 



b u = 2e" 



brr — 2 J, 



z z=0, 



^22; " ^rr? 

= 2r 2 efy,*|«= 0j 



& xx = 2e I V z | z=0 , 



-2e-^.U =0 



(3.70) 



(3.71) 



(3.72) 
(3.73) 



(3.74) 



(3.75) 



Including the extra dimensions the vierbein is replaced 
by the sechsbein, 



J 



\z) 



y/~9tt 



= 0,0,0, 



0,0,0,0,0 ,e (r) 



0,-=, 0,0, 0,0 ,e M J 

V 9rr 



= 0,0, 



'9<p<p 



:, 0,0,0 , 



0,0 1 0,0 ) -=,0) ,e (v) J 



0,0,0,0,0,- 



/<7i/j/ 



one can write down the energy-momentum tensor (|2,24p 



as 

[T A B }= - ee (t) A e {t) B +p r e {r) A e {r) B 

+ P v e {tp) A e (ip) B +p z e {z) A e iz) B 

+ Pxe( x) A e {x} B + p y e {y) A e {y) B , (3.76) 

yielding the energy density and pressures as 

f,z , 



Vr = [T r r ] = 0, 
Pz = [T Z Z ] = 0. 



P. 1 



Py = 0. 



(3.78) 

(3.79) 
(3.80) 
(3.81) 



Sirf 2 



(3.77) 



A second step to construct the disk is to choose two 
solutions of Laplace equations l|3.6ip and Q3,62p for the 
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functions cf> and v. Usually, for thin disks in AD it is com- 
mon to solve the problem with Schwarzschild or Chazy- 
Curzon solutions [30j (see Sec. Ill B[) . The £>-dimensional 
Schwarzschild solution has isometry group RxO(D-l). 
To write it in Weyl form, D — 2 orthogonal commuting 
Killing vector fields are required. For the Schwarzschild 
solution, this occurs only for D = 4, 5. Hence only the 
four and five-dimensional Schwarzschild solutions can be 
written in Weyl form. For D > 5, the D-dimensional 
Schwarzschild solution is not a generalized Weyl solution. 
However, the geometry obtained by taking products of 
the D = 4 or D = 5 Schwarzschild solution with asymp- 
totically flat space are easily seen to be Weyl solutions. 

We take the Newtonian potential associated to a rod 
of constant density (4fj| . 



(r, z) = In 



gi + Rj - 2 
Ri +R 2 + 2' 



(3.82) 



where i?i = ^Jf 2 + (z - l) 2 e R 2 = ^Jf 2 + (z + l) 2 , and 
f is the radial coordinate normalized by mass, f = r/m. 
The length of the rod is L = 2m. We also use a Chazy- 
Curzon solution for the v function (f(r, z) = [r 2 + (\z\ + 
a) 2 ] -1 / 2 ), which is asymptotically flat, resulting that our 
disk is in fact a Weyl solution. Thus, along the rest 
of the present work, the 4D part is solved using the 
Schwarzschild solution, or the Chazy-Curzon, and the ex- 
tradimensional part by using the Chazy-Curzon solution. 
Both pictures produce very similar results. 



where C x and C v are integration constants. 

Eqs. I|4.84p and (|4,86j) form a system of equations for 
(f 2 and i 2 . From these equations we find the rotation 
curves Vc, 



9<p<p dtp 
gtt dt 



V c = / ypv _ / z 



9w f 

gu i 2 



reduce to 



V c 



1 F{r)u. r + GfrQfrr 
F(r)v r - G(r)(2/r + <^, r ) : 



(4.87) 



(4.88) 



C~e 



where F(r) = -C^e~ v + C 2 y c y and G(r) = 1 - C 2 e 

Note that when C x = C y = (no extra di 
we have Vc 



mensions J , we nave Vc — \J ~ <j>. r /(2/r + <j)^ r ) that is the 
known formula for circular orbits in 4D Weyl geometry, 
Eq. I|2.39p . Another important task about the derived 
rotation curves is that Eqs. <|4.85|) and (|4.86[) provide 
that the disk actually is living in a universe endowed with 
Universal Extra Dimensions (a model where the fields can 
access the extra dimensions; for details about the model 
The boundary conditions in manner to 



see 



24, 25, 2 



and C y come from calculations 



determine values for C. 
of the stability of the disk, a subject studied in what 
follows. 



THE STABILITY 



IV. ROTATION CURVES 

From Eq. (|3.59[1 we have the first integral of motion, 

-e-H 2 + f(r 2 + z 2 )+r 2 c' t >ip 2 + c u x 2 + c-»y 2 = 1, (4.83) 

where x A = dx A /ds. Assuming f = and z = (par- 
ticles with no radial motion and confined on 2 = 0), 
Eq.JUHHJ reads 

- c-H 2 + r 2 eV + c v x 2 + c- u y 2 = 1. (4.84) 

By Eqs. (|3.77p - (|3.81|) . one evidently can see that extra 
dimensions pressures do not contribute to total disk pres- 
sures. Immediately one is tempted to argue that this is 
the same to say that extra dimensions do not contribute 
to the density profiles. However, Eq. (|3.63|) shows that 
f(r, z) contains extra dimensional field components and 
therefore a priori the rates x and y must be considered as 
non-null incognites. On this fashion, the geodesic equa- 
tions on the disk reduce to 

e v x = C x , e-"y = C y , (4.85) 
(e-*) tr i 2 - (rV),^ 2 
= C 2 x (e-n ir + C*(e v \ r , (4.86) 



We need to answer the question: what are the values 
for C x and C v (the integration constants) that make the 
disk stable? A classical method, the Rayleigh criterion, 
and an advanced one, a pertubative method, are both 
used to constrain a possible range of stable values for C x 
and C y . The correct rotation curves must be constructed 
following the results from a such stability. 

A. Rayleigh criterion 

The stability of the circular orbits in the disk plane 
can be studied using an extension of the Rayleigh sta- 
bility criterion [i^]. This method is extremely good for 
Newtonian systems, and only for circular orbits. Given 
a such limitation, here it will be used only as a reference 
method. In Rayleigh criterion, we have stability when 
h^p > 0, where h is the specific angular momentum of a 
particle in the disk plane (h = g lpip 'p), 

1 rC V 2r 2 e^(l/r + ^ ir ) 

(5.89) 
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We find that for different values of C x and C y less than 
unity stability is reached when a > 1. For small values of 
a (a < 1 , highly relativistic disks) we have a small zone of 
instability, typically around r = 3m. In Fig. |4]we show 
the stability of the disk by the criterion presented above. 
The principal constraint derived from this method is to 
hold suitable values for the cut parameter a. 



Let x A be an equatorial circular geodesic in a stationary 
axisymmetric space-time (|3.59p . i.e., the worldline x A = 
(t, r = const, ip = const+ftt, z = 0, x = const, y = const). 
Substituting the four velocity l|5.9ip and demanding that 
9ab,z (but not gAB,zz) vanishes in the equatorial plane, 
the components of Eq. (|5.90p for horizontal oscillations 
read: 



B. Perturbative method 

The stability of circular orbits in the disk plane can be 
studied using a more acurate and sophisticated method. 
In what follows we compute a perturbative method 
where we assume that the disk particles are describing 
equatorial circular geodesies in stationary axisymmet- 
ric fields. The perturbation of the geodesic equation 



T A bcx"x 



B^C 
A 



is done performing the transforma- 
tion x n — > x A + A A - where A A = (St, Sr, Sip, 5z, Sx, Sy) 
are infinitesimal elements. Therefore, equations for the 
perturbations are, 



A A + 2T A BC x B A c + T a BC mA d x b x c = 0, (5.90) 

where T A bc are the Christoffel symbols and x A are 
proper time derivatives dx A /ds and can be written for 
a circular orbital motion as 



x A = (u 4 , 0,0,^,0,, eg, 



(5.91) 



where u'f2 = Vc- Note that (|5.00|) is equivalent to the 
usual deviation equation [4l||. Assuming only horizontal 
oscillations in the 4D part of the disk (Sz — 0), we get 



(St,Sr,0,0,5tp,0,0). 



(5.92) 



(St) + 2T t tr u t (Sr) = 0, (5.93) 
{Sr) + 2T r tt + 2r w u*fi(^) + [(r«, r (5.94) 

r0 2 )(O 2 



'yy,rCy]Sr — 0, 



(Sip) + 2T ip vr nu t (Sr) = 0. 



(5.95) 



Suppose that the solutions for St, Sr and S(p have a form 
of harmonic oscillations, ~ e lKs , with a common proper 
angular frequency K. The condition for solvability of 
Eqs. |533l) -(|535 l) is then 



-K 2 




-K 2 





T r AB r u A u B 2iKT\ v u t n I = 0, 



2iKT^ ipr u t n 



-K 2 



(5.96) 

where T r A B, r u A u B = (T r tt , r + Y\ v ^ 2 )(u t ) 2 + 
Y r xx,rC 2 + Y r yy ,rC 2 . From the non-trivial solution of 
this equation we derive the oscillation frequency with re- 
spect to infinity K = K/u 1 , referred in literature as the 
epicyclic frequency [12], as 

k2 = o Z A P ( A + rfi + 3(f>,r/r + 3$ + Q(r)}, 
2 + r - F(r) 

(5.97) 

where Z(r) = — e ^/f, the metric is that given by 
Eq. (|3.59p and P(r) and Q(r) are terms related to ex- 
tradimensional imprints: 



P(r) = M(r)[2 + <f> >r r] + F(r)rv >r , 



(5.98) 



Q(r) = P(r)[H(r)/r + 0.5£ r /, r // - <j? r /2 - 0, rr /2] - H(r)<f>, r - 2H(r)/r - N(r)/r, 



(5.99) 



H(r) 



v-<i> 



-V — <f) 



2r 



-Ci\v, r f,rlf - K - V,rr] + -^^[-V r f, r /f - vf r + V. 



(5.100) 



N(r) = ^(r)[-3zv-2r^ r JV-O.5r 2 2 ri / r +O.5r 2 ^ 

(5.101) 



and where M(r) = C 2 e~ v + C 2 e v and F(r) = —C 2 e~ v + oscillations in Weyl geometry. Graphics for the squared 
C 2 e v . When C x = C y = ^> P(r) = Q(r) = (no extra epicyclic frequency is showed in Fig. [5l The configuration 
dimensions), Eq . (|5.97)1 becomes the known formula for is stable only when k 2 > 0. Our results show that the 
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FIG. 4: Stability of disks by the modified Rayleigh criterion. In (a) we show that stable disk cut parameters a are obtained 
only for a > 1 (the full line), where ft 4^ > 0. (b) In the region of interest, stable disks occur for extradimensional parameters 
C x = 0.8 (dotted line), C x = 0.85 (dashed line) and C x = 0.9 (gray line). For C x > 0.95, < 0, and the disk becomes 
unstable (full line), (c) In the region of interest, stable disks occur for extradimensional parameters C y = 0.1 (dotted line), 
C y = 0.2 (dashed line) and C y = 0.3 (gray line). For C y > 0.5, ft 4^ < 0, and the disk becomes unstable (full line). Here a is 
normalized by mass and r' = r/m. 
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FIG. 5: (a) Epicyclic frequency versus the normalized radius r' = r/m for a configuration where C x = 0.1 and C y = 0.85. The 
stability is achieved for < r' < 5.5, when n 2 > 0; a great number of rotation curves of spiral galaxies remains exactly into 
this region. The bold curve is the result for a 4D thin disk configuration, (b) Epicyclic frequency versus the normalized radius 
for a configuration where C x =0.1 and C y = 0.9. Here the stability is achieved for < r' < 15, when k 2 > 0. Here the disk 
cut parameter a is normalized by mass. We used r — r/m and a = 1.5. 
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TABLE I: Stable values for C x and C v 



Values for C x and C y 


Region where the disk is stable 


C — ft C — ft ("Npwtnnian'l 


All r' 


ft <r C <r ft 9 ft <r C ^04 


All r' 


ft < C <r ft 9 C — ft ^ 


< r' < 0.4 


< (7,. < 2 C„ — 7 


< r' < 1.4 


< C x < 0.2, C y = 0.75 


< r' < 2.6 


< C x < 0.2, C y = 0.8 


< r' < 4 


< C x < 0.2, C y = 0.85 


< r' < 5.5 


< C x < 0.2, C y = 0.9 


< r' < 15 


< C x < 0.2, C y > 0.95 


Unstable disk 



integration constants C x and C y have a very restricted 
range of stable values. In Table U it is possible to see the 
intervals for C x and C y where the disk is stable (here we 
fixed a = 1.5, a stable disk parameter according to the 
Rayleigh criterion; any values a > 1 are able to produce 
stable disks). 



A such stability study is fundamental to discuss what 
are the exact values to be used for the extradimensional 
parameters C x and C y . In general, extra dimensions con- 
tribute to destabilize the disk, but it is possible to derive 
a semi-phenomenological model where we have a very 
short range of stable values that can be used to fit ro- 
tation curves of galaxies or gravitational lensing of dark 
halos in a UED background (see Table [l| . Those sta- 
ble values can astrophysically constrain two-UED mod- 
els. A more detailed discussion on the present perturba- 
tive method is presented in 



C. Stable rotation curves 

The stability study of the previous section is funda- 
mental to discuss what are the exact values to be used 
for the extradimensional parameters C x and C y . In this 
sense, it is possible to derive a semi-phenomenological 
model where we have a very short range of stable values 
that can be used to compare to rotation curves of galax- 
ies. In Fig. [SK we show the rotation curves for some 
values of C x and C y , where we are not pondering yet the 
stable cases. In [6b we present only the stable curves, for 
various stable values of the cut parameter a. As discussed 
before, the stability is only achieved for < C x < 0.4 
and < C y < 0.95 (Rayleigh criterion) or the values 
from the perturbative method specified in Table [J (what 
restringes even more the range for C x and confirms the 



Rayleigh criterion for C y ). At the same time those val- 
ues also prevent super luminal behavior for the particles. 
When C x = C y = we have the usual AD general rela- 
tivists profile that is quite similar to a typical Newtonian 
one. On the other hand, the "flattening" of galactic rota- 
tion curves is an effect observed to occur in the galactic 
area that contains most of the baryonic (visible) matter. 
Beyond this area it becomes extremely difficult to ascer- 
tain the behavior of rotation profiles, so a model that 
just keeps these rotation profiles flat may be incorrect at 
larger cosmological scales. In Fig. we show that the 
curves derived are asymptotically flat and in cosmologi- 
cal scales the curves tend to zero. The density profiles as 
discussed previously are practically the same as in the AD 
case and we find that the extra dimensions affect smally 
the density and the azimuthal and radial pressures. The 
AD profiles can be seen in Fig. [2] and in Bicak et al. (3o| . 
In Fig. [7] we show a comparison between the profiles de- 
rived with and without extradimensional imprints, show- 
ing that the disagreement is small and the profile derived 
from the present model is similar to a AD one. Neverthe- 
less, it is exactly due to a such tiny disagreement, induced 
by the geometry of extra dimensions at large distances, 
that we have a correct profile for rotation curves. In 
other words, the presence of extra dimensions is affecting 
gravity at large distances. Maybe this effect is connected 
to theories on modified gravity [2l||, and the real nature 
of such modifications is due to extra dimensions. 

Taking into account only the stable curves, in Fig. [8] 
we compare with some optically observed rotation curves 
of spiral galaxies. Here we are not doing a composition 
of a halo dark matter velocities plus the disk gas and 
the velocities of stars. What is happening is that the 
clean stable general relativistic disk geodesies are sim- 
ply fitting the region of interest. The not surprising ad 
hoc adjustment of C x and C y actually could tell nothing 
about the astrophysical role of the extra dimensions in 
the model. However the calculation of stable disks brings 
over with it realizable values for C x and C y , what makes 
possible to visualize a minimum representation of a real 
disk galaxy. Those values produce the full line curves of 
Fig. [fib, fitting with great precision the region of interest 
(the plateau anomaly after r' ~ 2). We also compare to 
some phenomenological models used in astrophysics as 
Navarro, Frenk & White Q and the Courteau fit [H]. 
The observed data is taken from 44, 45, Hfj], but for 
an alternative data source see (4^1, where our model is 
also successful. Another important statement is that the 
present model recovers the Newtonian profile for r' — > oo 
(where the asymptotic function v — > 0, and the Newto- 
nian limit is reached). 

It is important to remember that this is not a pure 
phenomenological model, where we would be interested 
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FIG. 6: (a) Disk rotation curves with extra dimensional parameters C x = 0, C y = 0, i.e., usual standard relativistic profile 
(dotted line, stable), ; C x = 0.1, C y = 0.7 (dashed line, unstable); C x = 0.1, C y = 0.85 (gray line, stable); C x = 0.1, C y = 0.9 
(full line, stable), (b) Only the stable curves, for various values of the cut parameter a; here a is varying from a = 1 to a — 2.5, 
and C x — 0.1 and C y = 0.85. The dotted ones are the Newtonian-like (where there are no extra dimensions), and the full ones 
are the stable curves derived admitting two extra dimensions. In (c) we show that the curves derived are asymptotically flat, 
and in large scale they tend to zero. For all, the solutions used are Schwarzschild for the 4D part and Chazy-Curzon for extra 
dimensions. We take r' = r/m. 



in take a complete fit of observational curve. Here the 
interest is to achieve a reasonable explanation for the 
plateau anomaly using some of the stable calculated pa- 
rameters. In this sense, the model could be considered 
as a semi-phenomenological model. 



VI. RESULTS FOR SOLAR SYSTEM SCALES 

By simple quantitative arguments, one can easily show 
that, at solar system scales, the model presents the cor- 
rect observed results for rotation curves (a Newtonian- 
like behavior). The argument remains basically on the 
comparison between distance scales and the mass of the 
system. Consider, for simplicity, as an example based 
on observational real quantities, that a spiral galaxy has 
Mgai ~ lO 12 Af - 10 42 kg, with radius R ga i ~ 10 



kpc ~ 10 20 m. The normalization scale is then R' gal = 
RgailM ga i = 10" 22 m/kg. This gives, in Fig. Bt, that 
the r' scale (in natural units) for a such galaxy is around 
r' ~ 100 (where the observed rotation curve is around 
r' ~ 10). Now, doing the same estimative for solar 
system we have, considering roughly that the solar sys- 
tem mass is M so i ar ~ 10M Q ~ 10 31 kg and the radius 

R solar - 1 PC - 10 16 m, R' SQlar = Rsolar/ M so i ar - HT 15 

m/kg, or R' so[ar = 10 7 R' gal . Now, the equivalent r' ra- 
dius for solar system is r' = 10 9 . This last value remains 
in a absolutely Newtonian-like range for rotation curves 
(see Fig. Ejfc). Thus the immediate conclusion is that 
our model is in a correct fashion for both galaxies scales 
and the solar system scale, using the same stable parame- 
ters for C x and C y obtained from stability considerations. 
The model, or indeed extra dimensions, is then sensitive 
to R/M scales. 



14 








008 


>1 





006 


-H 






co 






a 





004 


Q 











002 











FIG. 7: Comparison of the density profiles of a disk consid- 
ering only AD dimensions (full line) and what is derived by 
the present example that considers 6D (dotted line). It is 
important to point that it is exactly due to this tiny disagree- 
ment that extra dimensions do affect rotation curves. Here 
a stable example with a = 1.5, for m = 1 (i.e. r' = r/m), 
C x =0.1 and C y = 0.85. Prom Eqs. §ZT7\ and (MB]) , the 
67? density depends on c/>,z/r', which clearly appears due to 
the influence of extra dimensions. For a large r', 4>, z /r' goes 
to zero and Eq. I|3.77p recovers the original AD profile. When 
r' — > 0, 4>, z /r' — > oo, but / goes faster to infinity, implicating 
that also Eq. I|3.77p recovers again the original AD profile for 
small r'. 



In Fig. [9] we show that all the curves of the model 
present the same result at solar system scale (the full line) 
and the dotted line is the Newtonian yl/V profile (or 
what is approximately observed at solar system scales). 
The difference between our model and what is expected 
for a Newtonian disk is relatively small (~ 10~ 5 ). 



VII. CLUSTER OF GALAXIES 

Here, using the simple example for 6D, we show that 
in the low acceleration regime, a universe endowed with 
UED predicts gravitational lensing of the correct magni- 
tude to explain the observations of intergalactic lensing. 
As in nature many elliptical galaxies and galaxy clusters 
are well modeled as spherically symmetric, we do our 



calculation for a spherically symmetric system plus two 
UED. We adopt the metric 



ds 2 



^dt 2 



e A [dR 2 + R 2 {d9 2 + sin 2 6>d^ 2 ) 
dx 2 +dy 2 }, (7.102) 



where, as before, x and y are the flat extradimensional 
coordinates. <& and A are functions only of R. Consider 
a light ray which propagates in the equatorial plane of 
the metric. The velocity x A of the ray must satisfy 



e*i 2 



-e A (ir + i?V+± 2 + ? /) =0. (7.103) 
From the stationarity of the metric follows the conser- 
vation law e*i = E where E is a constant characteris- 
tic of the ray. From spherical symmetry it follows that 
e A R 2 Lp = L where L is another constant property of the 
ray. Let us write R = (dR/dip)ip, x = (dx/dip)ip and 
y = (dy/dip)ip. Now eliminating t and ip from Eq. (|T. 103[) 
in favor of E and L, and dividing by E 2 yields 



(b/R)-e 
(dy/dtp) 2 



2 ~- A {R- 2 [{dR/dv) 2 



+ 1} = 0, 



{dx/dtp) 2 

(7.104) 



where b = L/E. By going to infinity where the metric 
factors approach unity one sees that b is just the impact 
parameter of the ray with respect to the matter AD dis- 
tribution center at R ~ 0. Rearranging the last equation 
we obtain the quadrature 



if = [1 + R 2 X + R 2 y ] 



211/2 



„A-4> 



1 -1/2 



~ -1 



d_R 
R ' 
(7.105) 

where R 2 = (dx/di?) 2 and R 2 = [dy/dR) 2 are the rates 
where the KK imprints are distributed along the AD 
galaxy cluster. At this point, where the physical met- 
ric exactly flat, this relation would describe a line with 
<f varying from to ir as R decreased from infinity to its 
value Rtum at the turning point, and then returned to 
infinity. Therefore the deflection of the ray due to gravity 
is 



A<p = [l + R 2 X 



R 2 y} 1/2 



Rtv 



-1/2 



dR 



(7.106) 



To shed some light in this last integral, one can benefit of the weakness of extragalactic fields which grant that A 
and $ are all small compared to unity. As consequence the above result is closely approximated by 

1/2 



(l + A-$) 



dR 



a = l 



(7.107) 
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FIG. 8: (a) Different models — our model (a Gravity with Extra Dimensions (GEDI) model) and [i^) for the rotation curves 
of the NGC 7331 spiral galaxy. Observational data from Sofue et al. (1999) [441 ] . (b) Different models (our model GEDI and 
for rotation curves of the high speed rotation spiral galaxy UGC 12591. Observational data from Rubin (1987) (45l |. (c) 
Different models (ou r GEDI model and [48Ll49l|) for the rotation curves of the NGC 3198 spiral galaxy. Observational data from 
Begeman (1999) [461 ] . The model fits with great precision the region of interest - the plateau anomaly after r/m ~ 3 in (a) 
and after r/m ~ 2 in (b) and (c). The stable parameters used are C x = 0.15 and C y = 0.88 for (a) and (b), and C x = 0.2 and 
C y — 0.85 for (c). For all, the solutions used are Schwarzschild for the 4D part and Chazy-Curzon for extra dimensions. The 
Navarro-Frenk- White [3] profile describes the halo curve and when combined with the disk gas curve it phenomenologically 
fits the rotation curve. The Courteau curve in (c) is only a reference of a pure phenomenological fit to rotation curves. 
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FIG. 9: Solar system scales. Comparison of rotation curves of a disk considering only AD dimensions, in a Newtonian y/T/r' 
solar system regime (dotted line, which coincides exactly to a AD disk regime — the dashed line) and what is derived by the 
present example that considers 6D (full line). Here a stable example with a = 1.5, r' = r/m, C x = 0.1 and C y = 0.85. For 
a 6D disk we have the same profile, but with greater velocities, although the difference is relatively small. The estimative for 
solar system, considering roughly that the solar system mass is M ao i ar ~ IOMq ~ 10 31 kg and the radius R so iar ~ 1 pc ~ 10 16 
m, R' ao i ar = Rsoiar /M so i ar = 10~ 15 m/kg, or equivalently, in natural units, r' radius for solar system is r' ~ 10 9 . 



The rewriting in terms of an a derivative allows us to Taylor expand the radical in the small quantity A — $ without 
incurring a divergence of the integral at its lower limit. The zeroth order of the expansion yields a well known integral 
which cancels the it. Thus, to first order in small quantities 



Aip ■■ 

And integrating by parts: 

2 [l+i?2+i?2]l/2 Q 



2 [l + i?2+i?2]l/2 q ,00 ( A _$) i?di? 



da J Wa (i?2_«62)i/2 



dc 



lim (A - $)(i? 2 

R— s-oo 



a=l 



(7.108) 



(7.109) 



Since $ and A decrease asymptotically as the in- 

tegrated term, being a independent, contributes noth- 
ing. Carrying out the a derivative, and introducing the 
usual Cartesian u coordinate along the initial ray by 
u = ±{R 2 -b 2 ) 1 / 2 , we have 

b [1 + R 2 r + R 2 } 1 / 2 f°° A n - $ „ 
A^=^ x - ' R - > R du. (7.110) 

A factor 1/2 appears because we have included the inte- 
gral in Eq. (|7. 109(1 twice, once with R decreasing to, and 
once with R increasing from b. The integral is now per- 
formed over an infinite straight line following the original 
ray. 

Fig. [10] shows that the new ray deflection using the 
above calculated model, for a range of projected rates 
R x and R y , is bigger than the deflection produced by 
a ray passing through a cluster calculated only with 4D 
general relativity (GR). The difference between a AD GR 



with dark matter and the above calculation (where the 
geometry of the problem substitutes the concept of DM) 
is that in a conventional dark matter scenario one would 
compute $ and A from Einstein equations including dark 
matter as source, whereas in the present model one has 
an aditional term and computes A and <& on the basis of 
the visible matter alone. 

Solving numerically the Einstein equations of the 
spherical metric and by a numerical integration of Eq. 
(|7. 1 10(1 . it is possible to obtain a comparison between a 
galactic cluster with and without extra dimensions. In 
Fig. [TTJwe demonstrate this, using the same values for C x 
and C y calculated in Sec. [V] In our model, for the cluster 
scales, extra dimensions geometry deflect more light as 
awaited. Thus the system is equivalent to a DM system 
where (79.5 ± 3.3)% of the total matter in a spherically 
symmetric galactic cluster is "dark". This is very similar 
to what happens accordingly to observations. It is veri- 
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FIG. 10: The new ray deflection using the above calculated 
model, for a range of projected rates R x and R y , is bigger than 
the deflection produced by a ray passing through a cluster 
calculated only with 4D general relativity (GR). The presence 
of extra dimensions acts exactly as a dark halo of cold dark 
matter. 

fied by lensing observations that cluster of galaxies 
are composed of three main components: ~ 5% in mass 
is the optically luminous baryonic matter in hundred of 
bright galaxies; ~ 10-15% is in the form of a bright X-ray 
inter-cluster gas; and the remaining ~ 80-85% is some 
sort of non-baryonic "missing mass". The explanation of 
our results is that such non-baryonic "missing mass" can 
be explained by the geometrical presence of Universal 
Extra Dimensions. One can picture the problem as the 
geometry provides an "equivalent" residual gravitational 
mass. The "extra" gravitational field is induced only by 
extra dimensions and observationally a such induction is 
seen as a missing mass. 

Actually, using simple quantitative arguments as in 
Sec. IVH one can see that a typical cluster of mass 
M clustcr ~ 10 15 A/ o ~ 10 45 kg and radius R c i usteI ~ 10 
Mpc - 10 23 m, has R/M scales of order R' ~ 10~ 22 
m/kg, similar to the values calculated for galaxies. In 
consequence, it is reasonable to expect a strange mod- 
ified behavior for cluster similar as those for galaxies. 
The final conclusion is that a multidimensional universe 
is sensitive to R/M scales. 

VIII. CONCLUDING REMARKS 

At the present paper, the possibility to construct a 
galactic disk embedded in a multidimensional space-time 
was investigated. Four principal results was achieved: i) 
The density profiles of the disk have a 4D similar shape, 
but present a very soft difference, indicating that extra 
dimensions are modifying the gravitational field at galaxy 
scales (r' ~ 1 — 100). ii) The rotation curves of stable 
disks fit the profile of observed galaxies, hi) The rotation 
curves for solar system scales are Newtonian-like, as ob- 
served astronomically, iv) A cluster of galaxies present 
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FIG. 11: Numerical integration of the gravitational lensing 
deflection produced by a cluster living in a universe endowed 
with UED. Here we use a equivalence between the gravita- 
tional field induced by extra dimensions and a fictitious mat- 
ter — "KK matter" — to show that the contribution of the 
"extra" field induced by extra dimensions acts as dark mat- 
ter for a galactic cluster. In the example above, the fictitious 
"KK matter" constitutes (79.5 ± 3.3)% of the total matter in 
a spherically symmetric galactic cluster. We used R y = 0.9 
and R x — 0.4, stable values according to the disk problem. 

a larger deflection for light rays than what is calculated 
for a AD cluster. 

To achieve such results, after demonstrate that vacuum 
Einstein equations can be written as a set of Laplace's 
equations for any number D of dimensions, we showed 
that if D is even, then the sum of extradimensional pres- 
sures cancel and the total pressure depends only of 4D 
components. Thus, the density profiles have practically 
the same shape as in the 4D case and we find that the 
extra dimensions affect little the azimuthal and radial 
pressures. As the simplest example, we calculated the 
case for a six dimensional space-time. Two constants of 
motion from projection of extradimensional particle ve- 
locities are the free parameters of the model. We had 
prevented the ad hoc adjustment of such parameters with 
observed rotation curves, preferring to investigate values 
where the disk becomes stable. The stability is achieved 
when the disk is Newtonian (where such parameters are 
null) or for a tiny range of values that makes the circular 
geodesies fit with great precision the rotation curves of 
many spiral galaxies. The stability calculation was done 
using both the Rayleigh criterion and a perturbative ap- 
proach. We compared such results to well succeeded as- 
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trophysical dark matter profiles and demonstrated that 
our predictions give the same gravitational lensing as 
does a dynamically successful dark halo model. We also 
introduced quantitative arguments to show that at solar 
system scales, the model presents the correct behavior 
for rotation curves (a Newtonian- like profile). 

The method to construct the disk is not new. Ac- 
tually we used well known general relativistic methods 
to construct thin disks, including, however, as a nov- 
elty, extra dimensions in the geometry of the system. 
Though the presented results introduce important con- 
straints to Kaluza-Klein theories. If it is true, maybe it 
would mean that the world we perceive is only a fraction 
of a greater invisible reality and maybe the 4D matter is 
a kind of modulation demanded both by the curvature of 
AD space-time and the extradimensional geometry. 

Two main questions could be formulated about the re- 
sults of the present paper: i) Do such results bring an 
alternative to the dark matter conjecture (showing that 
the strenght of gravitation may be modified by extra di- 
mensions, which would have similar consequences)? ii) 
Or do they rather suggest a particular form of dark mat- 
ter (consisting in Kaluza-Klein particles that may result 
from extra dimensions)? A first observation on this as- 
pect is that modifications in rotation curves come from 
the possibility that extra dimensions are endowed with 
dynamical test particles, i.e., the extra dimensions must 
contain the presence of fields (as similar to Universal Ex- 
tra Dimensions models). Eqs. H4.85J1 an d Q4.86)) show 
that the projection of extra dimensional rates is the de- 
terminant factor to obtain a deviation from a Newtonian 
profile. Becomes clear, from Sees. IVII and IVII[ that the 
effect of extra dimensions is to modify gravity according 
to the astronomical scale R/M (distance per mass), thus 
altering the rotation curves and lensing, for instance, in 
a more rigid framework than that of CDM models. In 
this sense, our model is similar to e.g. MOND or TeVeS 
2l| . but proposes and demonstrates that modifications 
of gravity comes naturally from the geometry of a multi- 
dimensional universe and that a multidimensional space- 
time is sensitive to R/M (distance per mass) scales. 
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Appendix: An alternative for the extradimensional 
modified gravity interpretation 



Alternatively, we can also present another interpreta- 
tion for the nature of the results. One can assume that 
actually the model is constraining a Kaluza-Klein dark 
matter particle (provenient from UED imprints) . See de- 
tailed discussions about a such topic in [25||. And maybe 
both interpretations (the pure geometrical approach and 
the induced KK matter), are physically equivalent (see 
for instance a idea proposed by Clifford in where 
not only could one describe fields through geometry, 
but also the particles that interact by means of those 
fields). About the particle aspect, although the main 
purpose of the present work is to construct a top-down 
model and not to shed light on a fundamental theory we 
can illustrate considering the possibility that our model 
could constrain a Kaluza-Klein dark matter particle to 
be tested at Large Hadron Collider (LHC) in next years. 
As an example, consider for simplicity that Kaluza-Klein 
dark matter particle in Universal Extra Dimensions could 
be considered as a fundamental particle, where is in- 
cluded also the possibility that a such particle could de- 
cay in others. Following |26j], one can derive a Lagrangian 
from the compactification of the extra dimensions and 
connect such compactification to the cross section of the 
particle. Using the generic notation x a , a = 0, 1, 3 + 6 
for the coordinates of the (4 + (5)-dimensional space-time 
(different from the usual non-compact space-time coor- 
dinates, x^, [i = 0,1,2,3), and the coordinates of the 
extra dimensions, y a , a = 1,...,S, the 4-dimensional La- 
grangian can be obtained by dimensional reduction from 
the (4 + <5)-dimensional theory, 
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C(xn = J d & y^-J^^Tz[Ff(x»,y a )F ia p{x»,y a )\ + C mges (x^y a ) 

+i (3,27,25) (x»,y a ) (T^D^ + T 3+a D 3+a ) (Q,U, V) T (a", y a ) 

Q(x»,y a ) {XuU{x^y a )ia 2 H*{x^y a ) + X v V{x^,y a )H{x^y a )) + h.c] } 



5.111) 



Here are the (4+(5)-dimensional gauge field strengths 
associated with the SU(3)c x SU(2)w X f/(l)y group, 
while £> M = d/dx^ - and Aj+a = d/dy a - -4 3+a are 
the covariant derivatives, with A a = — i Y^=i Si^ai^l 
being the (4 + <5)-dimensional gauge fields. The piece 
£mggs of the (4 + (5)-dimensional Lagrangian contains 
the kinetic term for the (4 + (5)-dimensional Higgs 
doublet H, and the Higgs potential. The (4 + 8)~ 
dimensional gauge couplings §i, and the Yukawa cou- 
plings collected in the 3x3 matrices Xu,t>, have dimension 
(mass) _5/2 . The fields Q,U and V describe the (4 + 5)~ 
dimensional fermions whose zero-modes are given by the 
4-dimensional standard model quarks. A summation over 
a generational index is implicit in Eq. I|8.11ip . For exam- 
ple, the 4-dimensional, third generation quarks may be 
written as Q< 0) = (t,b) L , uf ] = t R , vf ] = b R . The ki- 
netic and Yukawa terms for the weak-doublet and -singlet 
leptons, C and £, are not shown for brevity. The gamma 
matrices in (4 + 6) dimensions, r Q , are anti-commuting 
2 k+2 x 2 k+2 matrices, where k is an integer such that 
5 = 2k or 6 = 2k + 1. Chiral fermions exist only when 
5 is even, and correspond to the eigenvalues ±1 of T 4+s . 



The space-time is described as M 4 x T 5 , where the extra 
dimensions can be compactified in a T s tori. The 6D 
example of Section [TlIBI has a tori of radii r x = e v l 2 and 
r y = e~ v l 2 . Asymptotically the compactification tori ra- 
dius is constant and it is obtained by r c = ^Jr 2 + r 2 , and 

the periodic dimensions x k = x k + 2nr c , 6 — 6 < k < 5. 
The compactification is obtained by imposing the identi- 
fication of two pairs of adjacent sides of the torus retan- 
gle. In a pure UED model the compactication scale arises 
simply from the need to produce the right amount of DM. 
Here, a modified UED model, the compactification scale 
arises from the dynamics of galaxies and from clusters 
lensing. A detailed investigation about the compactifica- 
tion of the model will be done in future works, although 
our preliminary results indeed strength the possibility 
that our model could constrain, as one of the physical 
interpretations, a Kaluza-Klein dark matter particle to 
be tested at Large Hadron Collider (LHC), where the 
proposed energy range is 1-14 TeV, and the main signals 
of KK DM will be several tt resonances. 
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